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INTRODUCTION 


This  report  describes  activities  under  Contract  N00014-67-A-0 103-0007 
during  the  period  1 October  1973  to  Septeafcer  1974,  at  the  Department  of 
Atmospheric  Sciences,  University  of  Washington. 

During  the  past  year,  our  research  efforts  have  been  concentrated  In 
three  general  areas:  (1)  desalination  and  entrainment  processes  In  growing 

sea  Ice,  (11)  optical  properties  of  sea  Ice  and  consequent  effects  on  down- 
welling  irradiance  In  the  Arctic  Ocean,  and  ( 1 1 i ) theoretical  modeling. 
Laboratory  observations  were  made  of  brine  drainage  from  sea  Ice,  lateral 
melting  of  Ice  walls.  Interaction  of  petroleum  with  sea  Ice,  and  the  movement 
of  large  brine  pockets  under  strong  temperature  gradients.  Scientific  results 
from  previous  laboratory  work  on  the  growth  of  Ice  stalactites,  the  evolution 
of  under-ice  melt  ponds,  and  the  desalination  of  young  sea  Ice  have  been 
written  up  and  are  now  either  In  print  or  In  press.  In  connection  with  the 
radiation  In  Ice  program,  analysis  of  the  Barrow  data  has  been  completed  and 
a paper  describing  the  spectral  distribution  of  light  beneath  first-year  sea 
ice  has  been  submitted  to  a journal  for  formal  publication;  a paper  describing 
the  spectrophotometers  and  their  operatlona1  characteristics  h»„  also  bet.' 
published.  A field  program  to  Investigate  the  interaction  of  solar  radiation 
with  perennial  sea  Ice  and  the  ocean  was  carried  out  at  T-3  during  the  simmer. 
Spectral  data  on  albedos,  attenuation  within  the  Ice  and  snow,  and  the 
irradiance  field  beneath  the  ice  were  obtained.  In  the  area  of  numerical 
modeling,  rapid  progress  has  been  made  In  developing  the  theoretical  framework 
needed  to  combine  dynamic  and  thermodynamic  effects  In  a single  ice  model. 
Formulation  of  the  Ice  thickness  distribution  theory  is  essentially  complete 
and  has  been  incorporated  Into  the  AIDOEX  ice  model.  A complete  descrlptl;  i 
of  this  model,  together  with  the  results  from  construed  motion  calculations , 


-2- 


appeared  is  a joint  paper  In  a recent  AIDJEX  Bulletin.  In  addition,  a 
theoretical  study  of  the  conditions  under  which  Internal  wave  drag  effects 
could  be  important  to  the  aotlon  of  arctic  pack  Ice  has  been  made  and  a 

I 

description  of  the  results  published  In  the  AIOJEX  Bulletin.  Preliminary 
observations  downstream  from  a moderate  pressure  ndge  keel  near  T-3  during 
September  verified  the  presence  of  standing  internal  waves  near  the  pycno- 
cllne.  A detailed  field  study  of  such  waves  Is  planned  for  summer  1975. 

DYNAMIC  AND  THERMODYNAMIC  MODELING 

Work  In  this  area  has  concentrated  almost  entirely  upon  generalizing 
and  Improving  our  theoretical  treatment  of  the  ice  thickness  distribution 
(fi)  and  its  relationship  to  the  dynamics  and  thermodynaml ~.s  of  the  arctic 
Ice  pack.  A fundamental  problem  In  calculating  changes  In  the  thickness 
distribution  lies  In  describing  the  redistribution  of  Ice  during  ridge 
building.  Our  initial  treatment  was  one-dimensional,  in  that  It  considered 
only  the  effects  of  convergence  or  divergence  and  assumed  that  only  the 
thinnest  Ice  In  a given  region  went  Into  ridging.  However,  strain  measure- 
ments from  the  1972  AIDOEX  Pilot  Study  Indicated  that  shear  was  as  large  cr 
larger  than  the  divergence,  and  examination  of  'erlal  photographs  suggested 
that  shear  ridges  were  a common  feature  in  the  ice  pack.  Last  fall  the 
model  was  expanded  to  take  into  account  the  effects  of  shear  as  well  as 
convergence.  One  consequence  of  shearing  motions  is  that  ridges  can  be 
formed  from  Ice  of  almost  any  thickness,  and  therefore  the  thinnest  ice 
first  assumption  cannot  be  used  In  a two-dimensional  treatment.  Hence,  for 
I a given  state  of  strain,  i:  Is  necessary  to  specify  the  range  of  thicknesses 

I and  the  relative  amounts  of  Ice  consumed  In  the  ridging  process  together 

I with  how  this  Ice  Is  redistributed  Into  pressure  ice.  Few  field  data  are 
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available  on  this  subject;  however,  conservation  of  ar«a  and  conservation 
of  mass  place  strong  constraints  on  the  redlstrloutlon  process.  A general 
theoretical  framework  has  been  developed  which  Incorporates  thes?  constraints 
and  defines  the  range  of  acceptable  functions  which  characterize  the  ridging 
process . 

Although  we  obtained  a number  of  analytical  and  numerical  solutions  to 
the  thickness  distribution  equation  last  fall,  these  solutions  were  not 
realistic,  rather  they  served  only  to  provide  insight  into  the  behavior  of 
the  governing  equation  and  the  sensitivity  of  G to  our  basic  assumptions. 

What  was  lacking  was  a way  to  relate  the  thickness  distribution  to  the  large- 
scale  strength  of  the  ice.  This  was  finally  accomplished  by  D.  A.  Rothrock 
from  the  AIDJEX  Modeling  Group  through  conservation  of  energy  considerations. 
In  coor»erat1on  with  members  of  the  AIDJEX  Modeling  Group,  realistic  calcula- 
tions of  the  thickness  distribution  were  then  begun  using:  (i)  the  elastic- 

plastic  material  properties  of  the  ice  postulated  by  M.  D.  Coon,  (il)  thermo- 
dynamic growth  rates  predicted  by  the  Maykut/Unterstel ner  model,  end 
(111)  specified  strains  comparable  in  magnitude  to  those  observed  during 
the  AIDJEX  pilot  studies.  A series  of  cases  were  run  using  different  initial 
assumptions  regarding  the  shape  of  G for  a variety  of  strain  histories.  In 
addition  to  showing  how  the  thickness  distribution  changed  with  time,  the 
results  also  provided  time-dependent  Information  on  yield  strength,  potential 
energy,  and  the  plastic  work  performed  on  the  ice.  The  model  was  also  driven 
using  reel  strain  data  gathered  during  the  1972  AIDJEX  Pilot  Study;  results 
from  this  calculation  yielded  estimates  of  G and  the  stress  components  during 
the  forty  day  period.  The  calculations  demonstrated  that,  because  it  is  1 s 
distribution  of  the  thin  ice  which  determines  the  overall  strength  of  the  Ice 
pack,  thermodynamics  plays  a central  role  in  determining  the  material 


properties  of  the  Ice  and  Its  dynamic  behavior.  Since  growth  rates  define 
the  characteristics  of  the  Ice  thickness  equation,  the  thermodynamics  also 
determine  the  relative  Importance  of  the  Initial  conditions.  For  example, 
in  the  situations  considered  above,  it  was  found  that  the  yield  strength  of 
ti,«  Ice  was  strongly  dependent  upon  the  initial  choice  of  G for  the  first 
2.;-30  days,  after  which  time  It  was  primarily  dependent  on  the  strain 
history.  This  suggests  that,  in  the  winter  at  least,  it  may  not  be  as 
critical  to  know  the  exact  Ice  thickness  distribution  as  we  had  previously 
thou^tt.  A comprehensive  description  of  the  thickness  distribution  work 
and  the  elastic-plastic  approach,  together  with  results  from  the  construed 
motion  studies,  was  prepared  jointly  with  r>Mnbers  of  the  AIDJEX  Modeling 
Group  and  published  In  a recent  AIDJEX  Bulletin  under  the  title,  “Modeling 
the  Pack  Ice  as  an  Elastic-Plastic  Material". 

Recently,  we  have  used  the  wxl  to  analyze  macroscale  (1000  kilometer) 
strain  data  taken  from  the  triangle  formed  by  ARI.IS  II,  T-3,  and  NP-10  during 
the  two  year  period  from  May  1962  to  May  1964.  The  velocity  field  ox  the  ice 
for  this  case  proved  to  be  quite  similar  to  that  observed  on  a much  smaller 
scale  (100  kilometers)  during  the  1972  AIDJEX  Pilot  Experiment;  the  motion 
was  generally  diverging  with  the  shear  component  of  the  strain  tensor  being 
comparable  in  magnitude  with  the  divergence.  This  suggests  that  results 
obtained  on  the  100  kilometer  scale  chnser  for  the  AIDJEX  Experiment  «*ould 
Indeed  be  representative  of  the  large-scale  behavior  of  the  ice  pack.  Using 
growth  rates  derived  from  climatological  heat  flux  averages  together  with 
the  model,  we  were  able  to  calculate  how  the  thickness  distribution  within 
the  T-3/ARLIS  II/flP-10  array  changed  during  the  two  year  period.  Predicted 
seasonal  changes  in  G conformed  to  our  intuition  and  the  final  form  of  G was 
In  reasonable  agreement  with  submarine  data.  The  only  obv’ous  problem  which 
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showed  up  as  a result  of  the  calculations  was  that  the  model  apparently 
does  not  create  enough  thick  (greater  than  8-10  meters)  pressure  Ice. 

We  plan  to  look  more  carefully  at  the  problem  of  dynamic  redistribution  of 
the  Ice  to  try  and  cone  up  with  a rational  way  of  describing  the  production 
of  very  thick  pressure  ridges.  The  1000  kilometer  strain  data  were  also 
used  to  evaluate  the  sensitivity  of  the  model  to  assumptions  regarding  the 
redistribution  process  and  the  Initial  conditions.  In  general,  we  found 
that  the  yield  strength  underwent  large  fluctuations  in  time  and  was  signifi- 
cantly affected  by  our  initial  assumptions.  The  shape  of  G,  on  the  other 
hand,  was  onl>  rlldly  affected  by  these  assumptions  and  It  appears  that  the 
strain  history  Is  the  dominant  factor  in  determining  how  ice  thickness  is 
distributed.  A paper  describing  the  formal  development  of  the  thickness 
distribution  concepts  and  these  recent  results  is  being  prepared  jointly 
with  A.  S.  Thorndike  and  D.  A.  Rothrock  for  submission  to  the  Journal  of 
Geophysical  Research. 

The  thickness  distribution  calculations  carried  out  to  date  have  relied 
on  rather  crude  data,  however,  the  Main  AIDJEX  Experiment  should  provide  us 
with  a detailed  data  set  of  high  quality.  We  hope  to  be  able  to  use  these 
data  for  the  final  test  of  the  G model.  While  the  thermodynamic  growth 
rates  used  In  the  previous  calculations  may  be  considered  typical  of  Ice 
in  the  Central  Arctic,  there  is  no  reason  to  expect  that  conditions  over 
any  specified  area  In  the  ice  pack  will  conform  to  climatological  averages. 
For  a proper  test  of  the  model  we  must  have  growth  rates  which  are  represent- 
ative of  the  AIDJEX  strain  array.  Since  it  Is  impractical  to  measure  growth 
rates  under  all  ice  thicknesses  during  the  entire  14  month  period,  »e  have 
devised  a combined  theoretical-experimental  program  to  obtain  the  needed 
information.  AIDJEX  personnel  will  periodically  make  thickness  measurements 
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wnder  six  different  Ice  categories , Including  leads  and  pressure  Ice. 

These  results  will  be  used  to  verify  growth  rates  predicted  by  the  thermc- 
dynaeric  Model  on  the  basis  of  actual  teMperature  and  heat  flux  data. 

Because  the  turbulent  heat  exchange  over  thin  Ice  Is  substantially  larger 
than  that  over  thick  Ice*  It  will  be  necessary  to  alter  the  upper  boundary 
condition  In  the  Ice  growth  model  so  that  the  fluxes  of  sensible  and  latent 
heat  appear  explicitly  as  function  of  wind  speed  and  the  gradients  of 
temperature  and  moisture  in  the  lower  part  of  the  atmospheric  boundary  layer. 
Suitable  parameterizations  have  been  found,  but  we  have  not  yet  tested  them. 

One  problem  limiting  the  general  applicability  of  the  thermodynamic 
model  is  that  it  must  be  driven  by  sophisticated  heat  flux  data  which  are 
not  routinely  available.  For  thicker  ice  where  growth  rates  are  small  and 
not  very  sensitive  to  short  term  departures  from  the  mean,  this  Is  not  a 
serious  problem  since  climatological  averages  can  be  used  without  introducing 
large  errors.  This  is  not  the  case  for  thin  Ice  growth  and  it  would  be  useful 
tu  have  a model  which  utilizes  simpler  and  more  easily  obtainable  data.  To 
remedy  this  we  have  recently  reformulated  the  Ice  growth  model  so  that  it 
can  be  operated  using  only  data  on  air  temperature,  wind  speed,  and  snow  depth. 
This  "simplified"  model  would  be  used  primarily  to  describe  the  growth  of 
first  year  ice  and  would  not  be  capable  of  predicting  sunnier  ablation  at 
the  upper  surface.  It  would,  however,  represent  a great  improvement  over 
empirical  methods  of  predicting  ice  growth.  In  the  future  we  hope  to  be  able 
to  use  the  complete  model  to  parameterize,  to  some  extent,  the  transport  of 
heat  within  the  ice.  This  would  reduce  the  complexity  of  the  "simplified" 
model  still  further  so  that  It  could  be  coupled  directly  with  other  models 
of  the  ^ce  or  atmosphere. 
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EAOIATIOH  M ICE 

Analysis  of  the  Barrow  data  has  been  completed  and  the  transmission 
results  submitted  to  Limnology  and  Oceanography  under  the  title,  "The 
Spectral  Distribution  of  Light  beneath  First-Year  Sea  Ice  In  the  Arctic". 

This  paper  deals  primarily  with  the  amount  and  spectral  composition  of  the 
short-wave  radiation  which  would  typically  be  available  for  primary  production 
under  sea  Ice  In  the  coastal  portions  of  the  Arctic  Basin.  While  the  under- 
lying first-year  Ice  appears  to  be  fairly  homogeneous , surface  conditions 
were  found  to  have  a large  effect  on  total  transmission  through  Ice  of 
equivalent  thicknesses.  Three  main  types  of  Ice  were  considered  In  this 
study:  (1)  white  ice  which  was  above  the  local  water  table  and  was  continu- 

ously drained,  (11)  ice  covered  by  shallow  melt  ponds,  and  (Hi)  snow-covered 
ice.  The  radiative  properties  of  a fourth  type  of  ice,  blue  ice  which  was 
saturated  by  melt  water  at  the  local  water  table,  were  found  to  be  indistin- 
guishable from  those  of  the  melt  pond  Ice.  Examples  are  given  in  the  text 
showing  absolute  transmission  and  percent  transmission  as  a function  of  wave- 
length for  each  of  these  surface  types. 

Analysis  of  the  transmission  data  has  definitely  established  the  presence 
of  algae  within  the  lower  part  of  the  Ice.  Examination  of  the  data  sequence 
revealed  time-dependent  changes  In  the  shape  of  the  transmission  spectra 
which  we  suspected  were  due  to  decreased  algal  loading  as  the  underside  of 
the  ice  melted.  By  differencing  transmission  spectra  obtained  on  different 
days  beneath  Ice  of  similar  type  and  thickness,  we  were  able  to  determine  a 
measure  of  the  spectral  absorption  of  the  contaminant.  We  found  a strong 
absorption  peak  at  435  ran  with  a secondary  peak  at  670  nm.  The  curve  has 
almost  exactly  the  same  shape  as  those  obtained  by  Lorentzen  (1972)  In  the 
laboratory  with  cultured  brown  algae.  There  seems  to  be  little  doubt,  then. 
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that  there  were  algae  within  the  lower  part  of  the  Ice,  presifltably  living 
In  brine  channels  between  Ice  crystals.  Following  two  weeks  of  Intense 
melting  at  the  surface*  the  average  Ice  thickness  had  decreased  to  about 
120  cm  and  effects  fro «r  the  algae  were  no  longer  discernible  In  the  trans- 
mission spectra.  It  is  not  clear  whether  melting  at  the  bottom  or  flushing 
by  surface  melt  water  contributed  to  the  demise  of  the  algae,  but,  in  any 
case,  they  do  not  appear  to  be  able  to  withstand  vigorous  periods  of  melting. 

Using  the  Barrow  data  in  conjunction  with  a two-layer  photometric  model 
based  on  the  theory  of  Dunkle  and  Bevans,  we  were  able  to  obtain  predictions 
of  light  transmission  beneath  first-year  white  ice  and  melt  pond/blue  ice 
as  ? function  of  ice  thickness  (0-300  cm)  and  wavelength  (400-800  nm) . By 
assuming  that  the  spectral  distribution  of  the  incident  radiation  at  the 
surface  is  similar  to  that  observed  at  Barrow,  we  also  produced  contour 
plots  which  will  allow  cruder  estates  based  only  on  knowledge  of  total 
global  short-wave  radiation.  For  the  snow-covered  ice,  we  showed  total 
transmission  as  a function  of  both  snow  thickness  and  ice  thickness  Insofar 
as  possible,  all  effects  of  the  algae  were  removed  from  these  estimates. 

Thus,  the  graphs  should  allow  biological  oceanographers  to  calculate  the 
energy  available  for  primary  production  at  the  bottom  of  the  ice,  but  they 
provide  no  information  as  to  whether  that  producticr  takes  place  in  the 
ocean  or  the  ice. 

Because  of  the  extensive  preparations  necessary  for  the  s :n*ner  field 
trip,  we  did  not  have  time  to  prepare  a report  on  spectral  albedos  and 
attenuation  coefficients  in  first-year  sea  Ice  before  the  field  party  left 
for  T-3.  At  the  present  time  our  plans  are  to  combine  these  data  with 
similar  results  from  multi-year  ice  to  produce  a comprehensive  description 
of  the  optical  properties  of  arctic  sea  ice. 
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SlftHER  FIELD  PROGRAM  ] 

The  field  program  was  carried  out  at  T-3  during  a four  iron  period 
spanning  the  melt  season.  The  work  closely  followed  the  plan  set  forth 

i 

in  our  1974  proposal.  Most  of  the  observations  were  made  on  sea  ice, 

about  one  to  two  miles  from  the  edge  of  the  island.  The  actual  site  was  \ 

on  the  opposite  side  of  the  island  from  the  camp,  irfhlle  choosing  the  i 

site  this  far  from  camp  compounded  the  logistic  difficulties,  there  were 

a number  of  Important  advantages;  easy  access  to  the  sea  Ice,  minimal  , 

contamination  of  the  test  area  with  oil  and  particulate  matter  from  the  ] 

■\ 

main  camp,  and  sufficient  separation  from  the  island  to  avoid  anomalous  ] 

ice  thickening  and  water  runoff  problems.  Details  of  the  individual 

projects  are  as  follows:  ■ 

’.  Albedo  Observations  - Spectral  albedos  were  determ  ned  using  the  sub- 
mersible spectroohotcmeter  with  an  opal  glass  diffuser  to  measure  Incident 
and  reflected  Irradiance  directly.  A pair  of  Kipp  -ad i ^meters  was  employed 
to  simultaneously  obtain  total  albedos.  Albedo  data  were  taken  systemati- 
cally throughout  the  melt  season  and  provide  comparative  records  from  many 
locations  over  the  time  interval  when  the  surface  of  the  ice  pack  was 
undergoing  drastic  changes.  Supportive  measurements  of  grain  size,  snow 
density,  ond  thU  .ness  of  the  surface  scattering  layer  were  made  tc  assist 
the  interpretation  of  observed  albedo  changes.  Data  are  still  being 
reduced  but  we  find,  for  example,  that  white  ice  in  June  had  an  albedo  of 
about  .85  at  500  nm,  dropping  down  to  about  .40  at  1000  nm;  the  correspond- 
ing Kipp  albedo  was  close  to  .60.  For  wmdpacked  snow  the  albedo  at  500  nm 
was  just  above .90,  dropping  to  about  .80  at  1000  nm. 
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2.  Underlet  Radiation  - The  spectral  distribution  of  downwei 1 ing  Irradlance 
beneath  perennial  sea  Ice  was  aeasureu  with  the  submersible  spectrophoto- 
meter at  about  fifty-five  different  sites.  Data  were  recorded  beneath  white 
Ice  and  adjacent  melt  ponds  over  a thickness  range  of  115  an  to  250  cm.  The 
angular  distribution  of  the  light  beneath  the  Ice  was  measured  in  a number 
of  Instances  by  varying  the  attitude  of  the  Instrument.  He  found  that  the 
variation  In  light  Intensity  was  less  than  5%  over  a 45°  tilt  range.  Indica- 
ting that  the  light  beneath  thicker  ice  {at  least)  is  nearly  isotropic  in 
the  downward  hemisphere.  He  do  not  yet  know  whether  algae  were  present,  in 
our  test  area  during  the  underice  measurements.  He  saw  no  Indication  of 
algae  In  the  cores,  but  we  expect  the  shape  of  the  spectral  Irradlance  curves 
to  be  the  ultimate  test  of  their  presence.  Preliminary  extinction  coefficients 
for  white  Ice  have  been  calculated  and  show  a minimum  value  of  approximately 
.C074  cm-1  at  500  am,  suggesting  the  absence  of  algae  in  observable  quantities. 


3.  Attenuation  within  Ice  and  Snow  - Because  of  internal  melting  and  refreez- 
ing over  several  annual  cycles,  we  expect  multi-year  Ice  to  have  attenuation 
coefficients  which  are  substantially  different  than  those  of  homogeneous 
first-year  ice.  Multi-year  ice  typically  exhibits  a pronounced  scattering 
layer  near  the  surface  which  is  not  present  In  first-year  ice;  penetrating 
short-wave  radiation  also  modifies  Its  Internal  structure  giving  in  a different 
distribution  of  brine  volume,  bubble  density,  and  salinity  with  depth. 

During  the  sumner  field  experiment  we  attempted  to  determine  how  these 
factors  affect  the  optical  properties  of  the  ice.  Vertical  profiles  of 
forward-scattered  and  backward-scattered  radiation  were  taken  in  the  upper 
meter  of  the  Ice  and  snow  with  the  profiling  spectrophotometer.  Hhile  back- 
scatter  measurements  could  be  made  simply  by  Inserting  the  fiber  optics  probe 
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down  Mall  vertical  bore  holes,  the  forwardscatter  measurements  Involved 
sotawhat  more  effort.  Trenches  were  dug  In  the  Ice  and  spring  snow  and 
horizontal  holes  bored  about  1 meter  Into  the  side  walls  at  various  depths. 

The  fiber  optics  probe  was  then  mounted  so  that  It  looked  straight  up  and 
Inserted  to  the  end  of  each  of  these  holes.  In  the  sea  ice,  observations 
were  made  beneath  both  the  natural  decomposed  surface  layer  and 
an  artificially  scraped  surface.  Variations  in  grain  size  and  crystal 
structure  with  depth  were  also  recorded.  From  these  observations  we  should 
be  able  to  determine  the  extent  to  which  the  highly  scattering  surface 
layer  modifies  the  optical  properties  of  the  ice  and,  when  used  In  conjunc- 
tion with  the  underice  data,  to  infer  attenuation  coefficients  deep  within 
the  Ice. 

4.  Ablation  Program  - As  part  of  the  ablation  program,  a photographic  survey 
of  the  areal  coverage  of  melt  ponds  was  carried  on  throughout  the  melt  season. 
He  had  also  hoped  to  be  able  to  conduct  systematic  observations  of  ablation 
on  the  top  and  bottom  of  a moderately  large  pressure  ridge.  A suitable 
ridge  was  finally  located  and  an  array  of  ablation  stakes  set  up  across  the 
sail,  but  before  the  thickness  holes  could  be  drilled  a system  of  leads 
opened  and  separated  it  from  the  camp.  Because  access  to  the  area  was 
hazardous  and  time-consuming  and  because  there  were  no  other  suitable  ridges 
nearby,  It  was  decided  to  abandon  the  program.  This  data  Is  important, 
however,  to  large-scale  studies  of  the  mass  balance  of  the  Ice  pack  and  we 
have  therefore  written  it  up  as  part  of  a suggested  program  of  growth  rate 
observations  to  be  cerrled  out  during  AIDJEX. 
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5.  Internal  Ice  Structure  - The  objective  of  this  stu<j|y  Mas  to  monitor 
the  changes  In  Internal  Ice  structure  (brine  volume,  bubble  density,  and 
grain  size)  which  occur  as  a result  of  Interaction  with  solar  radiation. 
During  early  July  the  Ice  became  saturated  with  water  resulting  In  the 
development  of  a local  water  table,  usually  10-15  cm  below  the  surfa  e. 

Cores  from  melt  ponds  and  adjacent  bare  Ice  showed  no  appreciable  differ- 
ences In  either  brine  volume  or  bubble  density  below  the  water  table  during 
the  early  part  of  the  melt  season.  Above  the  water  table  the  Ice  was  white, 
granular,  and  partially  decomposed;  below  the  water  table  the  structure  of 
the  Ice  was  complex  with  Irregular  bubble  layers  appearing  throughout  the 
rores.  During  the  later  part  of  the  melt  season,  these  layers  had  a larger 
brine  volume  than  the  surrounding  Ice,  Indicating  selective  absorption  of 
short-wave  radiation.  Many  of  the  melt  ponds  exhibited  near-surface  bubble 
layers  which  ranged  In  thickness  from  1-15  cm  during  the  early  part  of  the 
melt  season.  All  but  the  thickest  of  these  bubble  sheets  were  ablated  away 
by  the  end  of  July.  Quantitative  measurements  Involving  slabbing  and  photo- 
graphing Ice  cores  began  in  August  when  the  samples  could  be  drained, 
prepared,  and  measured  without  appreciable  melting  in  the  open  air.  Cross- 
sections  1 to  2 cm  In  thickness  were  photographed  against  a background 
grid  of  known  dimensions  to  determine  bubble  densities  and  distributions  in 
the  upper  meter  of  the  ice;  samples  were  taken  from  beneath  melt  ponds, 
white  ice  and  newly  frozen  lead  Ice.  For  the  white  ice,  no  strong  systematic 
trend  with  depth  was  noted  below  the  local  water  table.  Values  ranged  from 
2-6%.  In  melt  pond  Ice,  however,  the  brine  volume  was  strongly  dependent  on 
depth  by  the  end  of  July;  at  the  end  of  the  melt  season,  the  liquid  volume 
ranged  from  40%  or  more  at  the  surface  down  to  10  or  15%  at  1 m.  Thus  the 
brine  volume  beneath  melt  ponds  Is  5 to  10  times  that  beneath  the  bare  Ice, 
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indlcatlng  a significant  storage  of  short-wave  energy  beneath  the  melt 
ponds  In  the  form  of  latent  heat. 

Spectrophotometer  data  were  routinely  recorded  both  on  magnetic  tape 
and  In  the  form  of  X-Y  plots.  The  X-Y  plots  were  Intended  to  function 
primarily  as  a backup  In  case  of  failure  In  the  magnetic  tape  records. 

Data  tapes  produced  In  the  early  part  of  the  melt  season  were  of  good 
quality  and  have  already  been  reduced t however,  the  tape  recorder  appears 
to  have  undergone  progressive  deterioration  during  the  course  of  the  sunnier 
and  later  tapes  show  an  Increasing  number  of  sporadic  coding  errors  and 
blank  spots.  By  the  end  of  the  experiment  the  tape  recorder  was  generating 
essentially  only  unrecognizable  characters  and  blanks.  Much  of  the  data 
from  the  middle  tapes  can  be  recovered  by  hand,  but  the  final  few  tapes 
contain  no  useful  Information.  Although  X-Y  plots  could  not  be  made  at 
every  site  due  to  occasional  plotter  malfunctions,  we  do  have  plots  from 
many  of  the  sites  where  the  magnetic  record  was  lost.  He  are  presently 
transcribing  magnetic  tape  data  from  the  middle  of  the  melt  season  for 
subsequent  keypunching  and  computer  reduction.  Because  of  these  processing 
problems  the  data  reduction  will  be  slower  than  originally  planned,  but  we 
should  ultimately  be  able  to  recover  a substantial  portion  of  the  missing 
data.  For  future  experimental  work  we  are  constructing  a state-of-the-art 
data  acquisition  unit  In  which  the  Information  Is  recorded  on  a cheap, 
battery  powered,  cassette  tape  recorder. 

The  radiation  related  observations  were  terminated  at  the  end  of  Angust 
and  the  month  of  September  devoted  to  a study  of  the  generation  of  internal 
waves  In  the  ocean  by  pressure  ridge  keels.  A theoretical  study  of  the 
problem  (see  next  section)  had  indicated  that,  under  the  proper  conditions, 
standing  Internal  waves  should  be  present  near  the  pycnocllne  downstream 
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from  pressure  ridge  keels.  No  previous  Measurements  of  these  waves  are 
known  to  have  been  Made.  At  the  beginning  of  the  study  we  were  not  opti- 
mistic about  our  chances  for  success  for  a number  of  reasons:  (1)  the 

part  of  the  ocean  where  T-3  was  located  was  fairly  stagnant  with  low 
relative  current  speeds,  (11)  the  only  large  ridge  (about  15  meters)  in 
the  area  was  the  one  chosen  for  the  ablation  measurements  which  was  still 
cut  off  from  the  Island,  necessitating  the  transport  of  equipment  across 
the  Ice  choked  leads  In  a small  boat,  and  (111)  we  would  be  able  to  get 
only  a rough  Idea  of  the  velocity  and  density  structure  of  the  ocean  near 
the  keel.  Nevertheless,  thermistors  were  lowered  to  the  pycnocllne  during 
the  first  week  of  September  and  observations  taken  for  a two  week  period. 

Internal  waves  were  detected  on  several  occasions,  with  the  largest  having 
an  amplitude  of  3-4  metert.  The  pycnocllne  was  relatively  shallow  1-  the 
area,  occurring  at  a depth  of  only  30  meters  as  compared  to  40-50  meters 
In  most  parts  of  the  Arctic  Ocean,  and  this  may  have  contributed  to  the 
unexpectedly  large  size  o*  the  waves.  While  the  data  were  too  crude  In 
time  and  space  to  determine  how  the  wavelength  and  amplitude  of  the  waves 
were  related  to  keel  geometry,  current  speed,  or  strength  of  the  pycnocllne, 
the  observations  suggest  that  such  waves  should  be  a frequent  occurrence 
throughout  the  Arctic  Basin.  Given  the  observed  ridge  frequency  of  several 
ridges  per  kilometer,  It  appears  likely  that  keel  generated  Internal  waves 
could  play  a significant  role  In  the  jpper  part  of  the  Arctic  Ocean,  and 
perhaps  In  the  motion  of  the  Ice  Itself. 

INTERNAL  WAVE  QUAG 

A simple  theoretical  model  describing  the  Internal  wave  drag  on  a 
pressure  ridge  keel  has  been  developed  by  F.  Rigby.  The  treatment  assumes 
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that  the  upper  pert  of  the  Arctic  Ocean  can  be  described  by  a two-layer 
model  with  a sharp  density  discontinuity  near  the  pycnocllne.  Mathemati- 
cally this  situation  Is  analogous  to  one  described  by  Land)  (1932)  Involving 
the  formation  of  waves  on  the  surface  of  a stream  which  are  induced  by  an 
obstacle  on  the  bed  of  the  stream.  With  some  modification  of  the  Laub 
approach,  anal,  leal  solutions  were  obtained  and  a description  of  the 
results  submitted  to  the  AIDJEX  Bulletin  under  the  title,  Theoretical 
Calculations  of  Internal  Wave  Drag  on  Sea  Ice". 

Calculations  were  carried  out  to  determine  the  dependence  of  Internal 
wave  drag  on:  (1)  depth  and  shape  of  the  keel,  (11)  relative  current 

velocity,  (111)  depth  of  the  Interface  between  the  two  layers  (h‘),  and 
(1v)  size  of  the  density  jump  across  the  Interface  (ftp).  Two  keel  shapes 
were  considered;  a semi-circular  keel  and  a semi -elliptical  keel  with  Its 
major  axis  three  times  as  long  as  Its  minor  axis.  This  semi -elliptical 
keel  was  very  similar  In  cross-section  to  one  studied  by  Kovacs  et  al.  (1973). 
The  results  show  that,  whereas  form  drag  on  the  semi -elliptical  keel  is 
substantially  smaller  than  on  the  semi-circular  keel,  wave  drag  at  high 
current  speeds  Is  much  larger  on  the  semi -elliptical  keel  than  on  the  semi- 
circular keel.  In  general,  the  results  indicated  that  wave  drag  is  neglig- 
ible below  current  speeds  of  about  25-30  cm/sec;  In  the  30-40  cm/sec  range, 
wave  drag  Increases  rapidly  and  can  become  the  dominant  component  of  the 
water  drag  for  the  larger  (>20  meters)  keels. 

Recently,  the  model  has  been  extended  to  treat  the  case  of  a homogeneous 
mixed  layer  of  constant  density  overlaying  a layer  with  a constant  density 
gradient.  This  Is  a closer  approximation  to  the  true  sltua^on  In  the  Arctic 
Ocean.  At  present  we  are  still  checking  over  the  details  o this  development 
and  have  not  compared  predictions  from  this  model  with  thotr  >f  the  simple 
two  density  model . 
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It  1$  evident  fro®  the  model  calculations  that  the  wave  drag  Is  very  sensi- 
tive to  assumptions  regarding  the  depth  of  the  Interface  and,  to  a lesser 
extent,  the  magnitude  of  the  density  jump  across  the  Interface.  At  this 
point.  It  Is  not  clear  what  values  for  h1  and  ip  are  most  appropriate. 

It  Is  Important  to  settle  this  question  as  the  applicability  of  our  results 
to  the  Arctic  depend  on  the  proper  choice  of  these  two  quantities;  however. 

It  appears  likely  that  wave  drag  Is  not  negligible  near  large  keels  when 
the  current  Is  flowing  strongly.  Even  If  wave  drag  proves  to  be  Important 
for  some  keels  In  some  situations.  It  Is  not  yet  clear  whether  this  effect 
would  appreciably  alter  large-scale  values  of  water  stress.  Direct  field 
observations  appear  to  be  the  only  practical  way  to  definitely  establish 
the  role  of  Internal  wave  drag  In  the  dynamic  behavior  of  the  Ice.  The 
preliminary  measurements  made  near  T-3  In  Septenfcer  1974  verified  the 
presence  of  these  waves  and  demonstrated  that  they  could  be  detected  with 
relatively  simple  and  Inexpensive  apparatus.  The  observations  also  suggested 
that  significant  Internal  waves  could  be  generated  without  extremely  high 
current  velocities.  Although  smaller  Internal  waves  might  not  greatly 
Influence  the  drag,  they  could  be  Important  to  understanding  processes  in 
the  upper  100  meters  of  the  Arctic  Ocean.  We  are  hoping  to  be  able  to 
conduct  detailed  observations  at  an  AIDJEX  camp  during  the  sumner  of  1975 
In  an  effort  to  obtain  more  quantitative  information  on  how  these  waves 
Interact  with  the  ocean  and  Ice  cover. 
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IABORATORV  STUOIES  Of  SEA  ICE 

During  this  reporting  period,  we  completed,  wrote  up,  and  edited 
the  results  of  the  following  four  projects.  First,  our  Ice  stalactite 
research  under  the  title  "Ice  stalactites:  Comparison  of  a laminar  flow 

theory  with  experiment"  by  S.  Martin  appeared  In  the  March  1974  Issue  of 
the  Journal  of  Fluid  Mechanics.  Second,  S.  Martin  and  P.  Kauffman  com- 
pleted their  joint  research  on  under- Ice  melt  ponds,  which  described  the 
observed  laboratory  formation  of  Ice  In  these  ponds  In  teres  of  three 
one-dlmenslonal  mathematical  models.  This  work,  under  the  title  "The 
evolution  of  under-ice  melt  ponds,  or  double  diffusion  at  the  freezing 
point",  appeared  In  the  July  1974  Issue  of  the  Journal  of  Fluid  Mechanics. 

Third,  S.  Martin  and  L.  Elde  analyzed  the  data  from  Elde's  experiments 
on  the  growth  and  desalination  of  a thin  sheet  of  laboratory  sea  Ice.  This 
material,  under  the  title  "The  formation  of  brine  drainage  features  In  young 
sea  Ice",  will  appear  shortly  In  the  Journal  of  Glaciology.  Finally, 

S.  Martin  and  W.  Campbell  wrote  a reply  to  a critique  of  their  earlier 
article  on  the  Interaction  of  oil  spills  with  the  large  scale  motion  of  the 
arctic  pack  Ice.  This  Is  scheduled  to  appear  In  the  29  November  1974  Issue 
of  Science. 

In  addition  to  the  work  described  In  the  above  papers,  our  present 
researci.  Is  concentrated  In  four  areas;  namely,  the  drainage  of  salt  from 
sea  ice,  the  lateral  melting  of  pressure  ridge  keels  and  other  ire  blocks, 
the  Interaction  of  petroleum  with  sea  Ice  and  the  migration  of  large  brine 
droplets.  He  will  next  discuss  these  four  in  the  above  order. 


1.  Drainage  Studies  - A third  year  graduate  student,  Mr.  Terren  Nledrauer, 
as  part  of  his  Ph.D.  research.  Is  doing  a laboratory  study  of  the  convective 
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drainage  of  brine  from  the  skeletal  layer  of  set  Ire.  His  apparatus  con- 
sists of  a small  freezing  cell  with  the  Inner  dimens  ons  of  10  cm  wide, 

20  an  deep,  and  only  1.5  m thi«k.  This  small  freezing  cell  Is  Insulated 
with  several  layers  of  thermopane,  while  the  flow  Inside  the  cell  Is 
recorded  by  time-lapse  movies,  still  photographs , and  small  nermocouples 
with  a thermal  t solution  of  10~2  degrees  Celsius. 

Using  this  apparatus,  Mr.  Nledrauer  observes  the  following  phenomena: 

(1)  A spatially  periodic  convective  drainage  pattern  which  occurs  between 
the  crystal  platelets.  The  pattern  consists  of  narrow  downward  jets  of 
brine  separated  by  broad  regions  of  upflow.  The  wavelength  of  the  pattern 
Is  about  5 m».  (2)  The  formation  of  brine  drainage  channels  from  the  join- 

ing together  of  large  brine  pockets.  (3)  The  expulsion  of  brine  from  very 
small  brine  tunes  with  a di  meter  of  hQ  microns  Into  larger  channels  and 
the  subsequent  closing  off  of  these  brine  tubes  Into  brine  pockets.  At  the 
present  time,  Mr.  Nledrauer  Is  both  working  on  theoretical  models  of  the 
above  phenomena  and  designing  as  a control  some  simple  three-dimensional 
experiments. 

2.  The  Ablation  of  Ice  In  the  Ocean  - The  problem  of  estimating  the  'lifetime' 
of  both  pressure  ridge  keels  and  Icebergs  In  the  ocean  involves  the  study  of 
the  laminar  and  turbulent  boundary  layers  which  occu*  on  both  vertical  and 
horizontal  Ice  walls  submerged  In  warm  sea  water.  During  the  past  six  months, 
a fourth  year  graduate  student,  Mr.  Edward  JostKugt.r,  who  1$  supported  by 
AIDJEX,  hes  been  studying  the  behavior  of  these  boundary  layers  both  theoret- 
ically and  In  the  laboratory  as  part  of  his  PhD.  research.  His  preliminary 
results  show  that  the  ablation  occurs  In  a complicated  double  convective 
boundary  layer. 
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To  describe  briefly  the  physics  of  this  double  boundary  layer,  con- 
sider a vertical  Ice  wall  In  a salt  water  solution  held  at  a temperature 
above  freezing.  At  the  Ice  wall,  the  boundary  condition  Is  that  the 
salinity  and  temperature  of  the  salt  water  Bust  lie  on  the  freezing  curve. 
Therefore,  the  wall  must  melt  to  cool  and  dilute  the  adjacent  warm  water 
down  to  the  freezing  curve.  From  a simple  steady-state  mathematical  model, 
because  heat  diffuses  much  faster  than  salt,  the  temperature  boundary  layer 
Is  about  3C  times  thicker  than  the  salinity  boundary  layer.  This  difference 
In  the  boundary  layer  thicknesses  has  the  following  effect,  depending  on 
whether  the  sea  water  salinity  Is  less  than  or  greater  than  about  25  °/oo. 

When  the  oceanic  salinity  Is  less  than  25  °/oo  and  the  oceanic  temperature 
is  below  the  temperature  of  maximum  density,  then  the  melting  wall  generates 
cold,  less-saline  water  relative  to  the  ocean.  The  less-saline  water 
generates  a narrow  upward  boundary  layer  jet,  while  the  colder,  thicker 
temperature  boundary  layer  also  yields  a broad  upward  jet,  so  that  the 
boundary  layer  flow  Is  one-directional. 

Idien  the  salinity  Is  greater  than  25  °/oo  and  the  sea  water  temperature 
Is  above  the  temperature  of  maximum  density,  then  the  salinity  boundary 
layer  still  flows  upward,  while  the  flow  In  the  temperature  boundary  layer 
is  now  downward,  so  that  a two-dl recti onal  flow  Is  produced  In  the  vertical. 
There  is  also  a gray  area  for  salinities  less  that  25  °/oo  and  temperatures 
greater  than  the  temperature  of  maximum  density  where  the  flow  regime  Is 
a combination  of  those  previously  discussed. 

Our  preliminary  experimental  results  show  for  the  first  case  described 
above  that  the  boundary  layer  becomes  turbulent  for  temperature  elevations 
above  the  freezing  point  of  order  1 degree  Celsius  In  a length  of  order  1 meter. 
The  transition  to  turbulence  Increases  the  melting  rate  by  about  a factor  of  5. 
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Durlng  the  cowing  year,  we  plan  to  study  both  the  laminar  and  turbulent 
regimes  for  the  various  cases  described  above  and  see  if  our  results  will 
lead  to  theoretical  Models.  Me  also  hope  to  apply  our  results  to  the 
proposed  observations  of  pressure  ridge  Melting  which  will  be  carried  out 
during  the  Main  AIOJEX  Experl went, 

3.  Oil  under  Sea  Ice  - Because  of  the  Increasing  possibility  that  off- 
shore oil  drilling  will  soon  begin  In  the  Beaufort  Sea,  and  the  subsequent 
Interaction  of  this  oil  with  the  sea  Ice  through  either  spills  or  blow-outs, 
we  have  done  several  slMple  experiments  on  how  sea  Ice  absorbs  oil.  Through 
the  offices  of  Imperial  Oil  In  Canada  and  Atlantlc-Rlehfleld  In  Alaska,  we 
obtained  free  samples  of  both  Norman  Wells  crude  oil  and  Prudhoe  Bay  crude 
oil.  Using  these  samples  and  diesel  fuel,  P.  Kauffran  and  S.  Martin  examined 
In  our  thin  Ice  tanks  how  the  entrainment  proceeds. 

Our  preliminary  experiments  show  that  a smell  fraction  of  the  oil  will 
be  entrained  by  the  brine  channels,  whereas  the  rest  of  the  oil  will  be 
trapped  In  a layer  at  the  Ice  front.  This  trapped  oil  apparently  has  several 
effects.  First,  because  the  thermal  conductivity  of  crude  oil  Is  about 
6 X 10*2  that  of  Ice,  tne  Ice  below  the  oil  layer  grows  more  slowly.  Second, 
the  oil  appears  to  completely  block  the  brine  drainage  from  the  sea  Ice  above 
the  oil  layer.  The  blocking  of  the  drainage  and  entrapment  of  brine  In  the 
Ice  above  the  oil  will  prevent  the  natural  Increase  of  both  strength  and 
conductivity  In  the  Ice.  Me  also  discovered  during  the  melting  of  these  Ice 
sheets  that  the  oil  rose  to  the  surface  through  the  brine  channels  which 
were  expanded  in  size  by  the  heat.  This  provides  a mechanism  for  bringing 
the  oil  to  the  Ice  surface  In  one  season.  P.  Kauffman  and  S.  Martin  are 
presently  attempting  to  organize  the  results  of  these  simple  experiments  for 
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publlcatlon.  Also,  S.  Martin  and  W.  Campbell  have  extended  their  comments 
on  the  Interaction  of  Urge  oil  spills  with  the  pack  Ice  dynamics  In  a 
recent  short  article  In  Science. 

4.  Migration  of  Large  Brine  Drops  - Experimental  studies  of  the  migration 
rates  of  macroscopic  brine  droplets  by  F.  Rigby  have  led  to  three  general 
conclusions,  first,  brine  droplets  with  a diameter  on  the  order  of  one 
millimeter  or  larger  can,  and  generally  do,  migrate  by  a convective  mecha- 
nism. The  droplets  display  speeds  as  much  as  two  orders  of  magnitude  greater 
than  those  predicted  theoretically  on  the  basis  of  salt  diffusion,  the  only 
mechanism  previously  considered.  Observed  velocities  ranged  from  0.2  to 
1.5  nn/hr.  Convective  brine  pocket  migration  thus  could  be  an  effective 
desalination  mechanism  In  sea  Ice  with  substantial  numbers  of  large  brine 
drops.  Secondly,  convective  brine  pocket  migration  is  characterized  by  a 
trail  of  flawed,  porous  ice  in  the  wake  of  the  moving  droplet.  The  trail  is 
physically  weak  and  represents  a path  along  which  subsequent  brine  droplets 
can  migrate  at  even  higher  rates.  This  observation  suggests  at  least  one 
possible  mechanism  by  which  brine  channels  can  Initially  form;  however,  there 
Is  no  experimental  evidence  to  Indicate  that  this  is  indeed  the  case.  Finally, 
the  data  Indicate  that  convection  takes  place  In  brine  pockets  only  if  the 
temperature  is  above  about  -5°C;  below  this  temperature,  diffusive  behavior 
is  observed.  Thus,  convective  brine  pocket  migration  can  be  expected  to 
occur  only  in  the  lower  part  of  the  ice  or  during  the  summer. 

Analysis  of  the  observational  data  has  been  hampered  by  the  fact  that 
the  results  do  not  appear  to  be  reproducible.  While  a particular  brine 
pocket  usually  migrates  in  a reasonably  consistent  fashion,  a second  drop 
placed  In  another  piece  of  Ice  under  similar  conditions  can  exhibit  < 
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rad  Icall  y different  rata  of  Migration.  It  appears  that  tha  Migration  Is 
aithar  vary  sensitive  to  snail  tanperature  fluctuations  or  to  Microfractures 
within  tha  block.  Careful  experlnents  are  needed  to  determine  the  origin 
of  the  velocity  variations,  but  it  seens  possible  that  the  density  of  cracks 
and  flaws  in  the  ice  nay  be  at  least  as  Important  as  the  site  of  the  tempera- 
ture gradient. 
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Recent  field  observations  in  the  polar  oceans  show  that  the  hollow 
tubes  of  Ice  called  ice  stalactites  form  around  streamer:  of  cold 
orine  rejected  by  the  growing  sea  ice.  In  a laboratory  study  of 
this  process,  we  inject  cold,  dense  brine  at  a constant  salinity, 
temperature  and  volume  flux  into  an  Insulated  tank  of  sea  water 
held  at  Its  freezing  point,  then  photograph  the  resultant  stalac- 
tite growth.  Because  the  inner  wall  temperature  of  the  stalactite 
remains  on  the  saiinlty-determlned  freezing  curve,  as  the  stalac- 
tite grows  and  the  temperature  deficit  of  the  brine  goes  into  the 
growth  of  Ice,  the  inner  wall  melts  to  dilute  and  cool  the  adjacent 
brine  back  to  its  freezing  point.  This  melting  means  that  both  the 
inner  and  outer  stalactite  radii  increase  with  time.  The  radius  of 
the  stalactite  tip,  which  is  constant  for  each  experiment,  is 
shown  to  be  controlled  by  the  onset  of  a convective  instability. 

If  fht  tip  become:  too  large,  overturning  occurs  and  the  sea-water 
intrusion  freezes,  reducing  the  radius  of  the  tip  so  that  the  flow 
leaving  the  tip  is  marginally  stable.  Inside  the  stalactite, 
since  the  inner  radius  increases  with  time,  both  theory  and  experi- 
ment show  the  interior  flow  to  be  convectively  unstable.  The 
present  study  also  derives  a solution  from  the  constant-heat-flux 
Graetz  solution  for  the  grows  In  both  length  and  side-wall  area 
of  the  stalactite.  The  experiments  show  that  away  from  the 
stalactite  base  c.id  the  very  beginning  of  the  experiment  this 
solution,  wi in  convection  accounted  for  by  an  adjustable  coefficient, 
describes  the  experimental  growth.  Finally,  analysis  of  the 
experiments  shows  that  as  much  as  50%  of  the  ice  represented  by 
the  cold  brine  does  not  go  into  the  stalactite,  rather  the  ice  goes 
directly  into  the  ocean  as  loose  crystals. 
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A model  of  the  motion  of  drifting  pavA  ice  in  the  Arctic  Ocean  is 
presented,  treating  explicitly,  If  not  with  consistent  rigcr,  the 
growth  and  melt  rates  of  the  Ice,  the  formation  of  leads  and  pressure 
ridges,  and  a mechanical  response  which  is  elastic  at  low  stress 
levels  ana  plastic  at  some  highe*\  critical  state  of  stress. 


The  strength  of  the  ice  Is  determined  by  Its  thickness  distribution, 
and  therefore  varies  because  of  both  thermal  and  mechanical  effects. 
To  examine  the  behavior  of  the  model,  several  artificial  calcula- 
tions were  made  by  specifying  the  strain  rate  history  of  a single 
element  of  pack  Ice  and  solving  for  the  Ice  thickness  distribution 
and  the  states  of  stress  In  the  Ice. 
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Strain  and  strain  rate  calculations  are  presented  based  on  measure-  t 

raents  of  the  positions  of  three  drifting  sea  Ice  stations.  The  j 

measurements  were  made  using  the  Navy  Navigation  Satellite  System.  ^ 

In  processing  the  data,  time  derivatives  were  calculated  by  Kalman  ? 

smoothing  techniques,  and  space  derivatives  by  simple  finite  ? 

differences.  ] 
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In  an  experimental  and  theoretical  study,  we  model  a phenomenon  * 
observed  in  the  summer  Arctic,  where  a fresh-water  layer  at  a | 
temperature  of  0°C  floats  both  over  a sea-water  layer  at  its  freezing  | 
point  and  under  an  ice  layer.  Our  results  show  that  the  1c  growth  ; 
In  this  system  takes  place  In  three  phases.  First,  because  the  fresh-  i 
water  density  decreases  upon  supercooling,  the  rapid  diffusion  of  § 


heat  relative  to  salt  from  the  fresh  to  the  salt  water  causes  a 
density  Inversion  and  thereby  generates  a high  Rayleigh  number 
convection  In  the  fresh  water.  In  this  convection,  supercooled 
water  rises  to  the  ice  layer,  where  It  nucleates  into  thin  vertical 
Interlocking  ice  crystals.  When  these  sheets  grow  down  to  the  inter- 
face, supercooling  ceases.  Second,  the  presence  of  the  vertical  ice 
sheets  both  constrains  the  temperature  T and  salinity  s to  lie  on 
the  freezing  curve  and  allows  them  to  diffuse  In  the  vertical.  In 
the  Interfacial  region,  the  combination  of  these  processes  generates 
a lateral  crystal  growth,  which  continues  until  a horizontal  ice 
sheet  forms.  Third,  because  of  the  T and  s gradients  in  the  sea 
water  below  this  Ice  sheet,  the  horizontal  sheet  both  migrates  upwards 
and  Increases  in  thickness.  From  one-dimensional  theoretical  models 
of  the  first  two  phases,  we  find  that  the  heat-transfer  rates  are 
5-10  times  these  calculated  for  classic  thermal  diffusion. 
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5.  Roulet,  R.  R.,  G.  A.  Maykut,  and  T.  C.  Grenfell,  Spectrophotometers  for 
the  measurement  of  light  In  polar  Ice  and  snow.  Journal  of  Applied 
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Two  portable  spectrophotometers  have  been  designed  to  record  light 
scattering  and  absorption  in  polar  Ice  and  snow.  In  the  first 
Instrument  optical  fibers  are  used  to  transmit  light  from  the 
Interior  of  the  Ice  to  the  spectrophotometer.  Such  an  arrangement 
allows  light  measurements  up  to  2 m away  from  the  instrumentation 
with  minimal  disturbance  of  the  natural  environment.  A miniaturized, 
submersible  spectrophotometer  was  also  built  for  In  situ  measurements 
under  floating  sea  Ice.  This  version,  except  for  the  recording 
apparatus.  Is  entirely  self-contained  and  is  housed  in  a cylindrical 
tube  9 cm  In  diameter  and  60  cm  In  length.  The  unit  can  be  lowered 
Into  the  ocean  through  a small  borehole  in  the  Ice;  position  and 
orientation  are  controlled  from  the  surface.  Both  spectrophoto- 
meters are  designed  to  measure  light  intensities  In  the  visible 
spectrum  (400-1000  nm).  Wavelength  resolution  is  adjustable  down 
tc  C nm  at  a wavelength  of  400  run,  with  a field  of  view  of  less 
than  3°.  Sensitivities  in  the  present  versions  are  sufficient  for 
measurements  through  several  meters  of  sea  Ice  with  a relative 
accuracy  of  IX.  Instrument  operation  has  been  tested  in  the  Arctic 
down  to  temperatures  of  -25°C. 
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Because  of  the  density  stratification  in  the  upper  part  of  the 
Arctic  Ocean,  the  movement  of  pressure  ridge  keels  through  the 
mixed  layer  can  create  <'*ernal  waves  in  the  vicinity  of  the 
pycnocline.  These  waves  cransport  energy  away  from  the  keels 
and  thereby  generate  drag  on  the  ice.  This  study  defines  the 
conditions  under  which  wave  drag  could  contribute  significantly 
to  the  total  water  stress.  The  development  of  Internal  waves 
downstream  from  a semi  elliptical  keel  is  described  by  a simple 
two-dimensional,  two-layer  ocean  model.  Results  calculated  from 
the  model  show  the  dependence  of  wave  drag  on  keel  depth,  current 
speed,  depth  of  the  upper  layer,  and  density  change  across  the 
interface  between  the  layers.  Comparison  of  the  relative  magni- 
tudes of  form  drag  and  wave  drag  suggests  that  wave  drag  is  not 
negligible  near  large  keels  when  the  current  is  flowing  strongly. 
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